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ABSTRACT. This paper outlines a method for reducing astrometric data to derive the closest approach 
time and distance to the center of an occultation shadow for a single observer. The method applies to CCD 
frames, strip scans or photographic plates and uses a set of field stars of unknown positions to define a 
common coordinate system for all frames. The motion of the occulting body is used to establish the 
transformation between this common coordinate system and the celestial coordinate system of the body s 
ephemeris. This method is demonstrated by application to the Tr60 occultation by Triton on 1993 July 10 
UT. Over an interval of four nights that included the occultation time, 80 frames of Triton and Tr60 were 
taken near the meridian with the U.S. Naval Observatory (USNO) 61 -inch astrometric reflector. Application 
of the method presented here to these data yields a closest approach distance of 359±133 km 
(corresponding to 0.017±0.006 arcsec) for the occultation chord obtained with the Kuiper Airborne 
Observatory (KAO). Comparison of the astrometric closest approach time with the KAO light-curve 
midtime shows a difference of 2.2±4.1 s. Relative photometry of Triton and Tr60, needed for photometric 
calibration of the occultation light curve, is also presented. 


1. INTRODUCTION 

From an occultation observation, one can learn about 
planetary rings and atmospheres with kilometer-scale spatial 
resolution (Elliot 1979), but to take full advantage of these 
data requires accurate knowledge of the location of the ob- 
served chord relative to the center of the occultation shadow. 
This must be known with an accuracy (ideally better than a 
few kilometers) such that the uncertainty introduced by as- 
trometric errors into the physical parameters derived from 
the light curve is smaller than their formal errors from model 
fitting. The most reliable method for determining the location 
of an occultation chord relative to the center of the shadow 
has been to use the immersion and emersion times from a set 
of occultation chords to fit a two-dimensional figure to the 
body shadow. For the giant planet atmospheres, which are 
rapidly rotating, an oblate figure is fit to the “half-light” 
times (French et al. 1985; Baron et al. 1989). For a slowly 
rotating atmosphere, such as Pluto, a spherical model is ad- 
equate (Millis et al. 1993). For large, airless bodies, such as 
Pallas, an elliptical figure has been used (Wasserman et al. 
1979), but with a larger number of chords, more details can 
be inferred about the figure of the body (Dunham et al. 
1990). 

Other methods for determining the location of the occul- 
tation chord relative to the center of the occultation shadow 
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include (for Saturn and Uranus) using ring occultations to 
establish the astrometry and learn where in the atmosphere 
the chord probed (Baron et al. 1989; Nicholson et al. 1995), 
since the orbits of the rings are known quite accurately 
(French et al. 1988; Elliot et al. 1993). Also, a central flash 
can be used to establish the position of a single chord 
(Lellouch et al. 1986). Finally, one can use the length of the 
chord along with the known figure of the body at a specific 
pressure level (Sicardy et al. 1991) to achieve the same end. 
However, in using that method, one must assume that the 
atmospheric structure is the same at the time for which the 
figure was established and at the time of the single-chord 
observation. 

For the Triton occultation presented here, we have only 
one observed chord (from the KAO), so the standard method 
of fitting immersion and emersion times from multiple 
chords cannot be used for the astrometric calibration. We 
were fortunate to obtain high-quality astrometric data to pre- 
dict the occultation. These data produced such an accurate 
prediction that we were led to investigate just how well they 
could be used in a post-event analysis, especially since addi- 
tional data were recorded the night after the occultation as 
well. 

We refer to the method presented here as the ephemeris 
method for reasons that will become apparent. We make the 
basic assumption that stars do not move significantly over a 
time span of a few days and use a common set of field stars 
to define a reference system (in detector coordinates) across 
all frames. We can use the occulting body’s ephemeris to set 
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the scale in right ascension and declination. With this ap- 
proach we do not have to reconcile two or more systems of 
celestial coordinates (the occulting body’s ephemeris and 
that of an astrometric network), and we are immune to scale 
errors that must be present at some level in an astrometric 
network. 

We shall first describe the method for determining the 
placement of the observed chord relative to the center of the 
occulting body from astrometric data. The next section de- 
scribes the observations and details of the astrometric and 
photometric calibration for the Tr60 occultation (McDonald 
and Elliot 1992). Finally, we present a discussion of the error 
and considerations for application of this method to future 
occultations. 

2. OVERVIEW OF THE EPHEMERIS METHOD 


the reference system. Even if some frames are missing ref- 
erence stars (due to poor seeing, or inconsistent centering 
causing stars at the edges to be lost), the mean frame can still 
be constructed. Each frame is linearly registered to the de- 
tector network to transform the occulted star and occulting 
body centers to a common reference system. The origin of 
this system is shifted to be the weighted mean center of the 
occulted star. Next, the scale and center of the detector net- 
work are related to celestial coordinates through a linear reg- 
istration of the occulting body centers to the ephemeris. This 
allows us to determine the right ascension and declination of 
the occulted star in the coordinate system defined by the 
body’s ephemeris. From here, we can determine the closest 
approach of the occulting body to the star and therefore, the 
placement of the occultation chord to the center of the body’s 
shadow. 


The method outlined here can be applied to CCD stare 
frames, CCD strip scans or photographic plates. Some of the 
frames (or plates) should contain the occulting body and oc- 
culted star both before and after the occultation event. Addi- 
tional frames, without the occulting body, can be used to 
reduce the error in the occulted star center, but this is not 
typically the leading source of error. Also, for best results 
(highest precision) the mean center of the occulting body 
observations should be as close to the occulted star center as 
possible. 

This paper deals solely with the reduction of stare frames. 
This analysis could be applied to strip scans, but the low 
frequency errors in the star positions (Dunham et al. 1991) 
would have to be removed first. Strip scan data used for the 
occultation prediction were not used in this analysis because 
of a field distortion found in the data (Dunham, personal 
communication). 

The ephemeris method is based on three assumptions. The 
first is that the field stars used for the common reference 
system are fixed (over the time scale of the observations — a 
few days). The second assumption is that the occulting 
body’s motion can be used to set the scale in right ascension 
and declination. The ephemeris does not have to be correct in 
terms of absolute coordinates because we are only interested 
in relative positions. The final assumption is that displace- 
ments of observed positions of the stars and occulting body 
relative to the mean centers (the centers of the stars without 
displacing effects such as refraction and aberration) can be 
removed by a linear registration. We will break up the 
sources of these displacements into two categories: (i) 
known effects such as refraction, diurnal aberration, annual 
aberration, and general relativistic bending, and (ii) unknown 
nonlinear effects such as field distortion of the telescope, 
detector distortion, nonuniform thermal distortion or signifi- 
cant chromatic differential refraction. 

In the following discussion, we use the term registration 
of a set of positions to a coordinate system to mean a least- 
squares fit that determines the linear transformation coeffi- 
cients [which we shall also call registration coefficients ; see 
Eq. (2) below]. 

This analysis is based on the construction of a common 
reference system— -the detector network, which is the 
weighted average of the centers of each field star chosen for 


3. DESCRIPTION OF THE METHOD 

In this section we present the method in detail. First, one 
establishes a set of stars that are common to all Nj frames 
and determines their centers by a suitable procedure. To be 
definite, we shall refer to CCD data and the row and column 
center of the nth star in the &th frame as (r nk ,c nk ). For 
photographic plates, one might choose a different notation. 
For the &th frame, the row and column coordinates of the 
occulted star are denoted by and those for the 

occulting body (in many cases, a “planet”) by (r pk ,c pk ). 


3.1 Creating the Detector Network 


A common coordinate system for all the frames is needed, 
so a detector network is created from the weighted average 
position for each of the selected field stars. We weight the 
data by frame, since the frames are not necessarily the same 
quality due to varying observing conditions (seeing, air- 
mass). For the data used here, the errors in the star positions 
appear to be the same for a single frame, independent of 
stellar magnitude (over the magnitude range of 12.6 to 18.7 
from a yellow photographic emulsion, kindly supplied by A. 
Klemola). For this case, the weights, w rk for the rows and 
w ck for the columns, are derived from a registration of each 
frame to an unweighted detector network. For the nth star we 
find its coordinates in the unweighted detector network by 
the weighted averages (with all the weights set equal to 1): 


^ k L^rk^nk 

2 *ii w v* 


S'Mf 

z k = | Viijt 

* 


(1) 


Then each frame is registered to the unweighted detector 
network, so that the reciprocal of the variance of the row and 
column residuals from the least-squares fit can be used as the 
row and column weights for each frame, and Eq. (1) is used 
a second time. The result is the detector network. 

For each frame, the row and column positions of the net- 
work stars are registered to the detector network to u^rmine 
the registration coefficients, a . 
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( 2 ) 


The observed row and column centers for the occulting 
body and occulted star are transformed to the detector- 
network coordinate system with the registration coefficients 
found in Eq. (2). The resulting occulting body coordinates 
for the fcth frame are denoted by (r rpk ,c rpk ). Similarly, the 
occulted star coordinates in this detector-network system are 
(r r * k ,C r * k ). 


3.2 Relating the Mean Frame to Celestial Coordinates 

Since the row-column coordinates of the detector are in a 
plane, we convert the celestial coordinates of right ascension 
and declination (RA, a, and Dec. S) to those in a plane 
tangent to the celestial sphere (assumed parallel to the detec- 
tor). The optical axis of the telescope intersects the mean- 
frame coordinate system at the center of the detector which 
correspond to celestial coordinates (o^.Sq). Here we intro- 
duce the tangent plane , as defined by Smart (1977), that has 
the coordinate £ parallel to RA and 77 parallel to Dec. The 
conversion equations between (a, 8) and (£ rj ) are (Smart 
1977): 


cos £sin(a— a 0 ) 

sin Sq sin £+cos Sq cos S cos(a- a 0 ) 5 

cos Sq sin S— sin 3$ cos £cos(a-a 0 ) 
sin Sq sin S+co$ cos £cos(a-a 0 )’ 


(3) 


As previously mentioned, we use the body’s ephemeris to 
establish the (a,S) system. The ephemeris should be topocen- 
tric and include a correction for light-travel time, but may or 
may not contain corrections for diurnal and annual aberration 
(see Sec. 4). Due to our third assumption (that the effects of 
refraction, diurnal aberration, annual aberration, and other 
displacements of the star centers can be removed by a linear 
registration), these corrections can either be included or not, 
since they will not make a difference in the resulting relative 
position of the occulting body and occulted star. 

We generate ( a eJc ,S ek ) from the ephemeris at the mid- 
times t k , of the frames where the subscript “e” stands for 
ephemeris. From ( a ek , S ek ) we use Eq. (3) to calculate the 
tangent-plane coordinates (£«,*, Ve*) f° r the occulting body 
as given by the ephemeris. The registration of these tangent- 
plane coordinates to the occulting body’s positions in the 
detector-network coordinate system provide the transforma- 
tion to celestial coordinates (as defined by the ephemeris of 
the body). 

We center this transformation on the occulted star. To do 
this we form the weighted average row and column position 
in the detector-network system. We denote this center by 
i(r r *k),(c r *k)), and it is computed with the weights 
(w rk ,w cA ). This shift of origin does not change the solution; 
it only simplifies later equations [Eqs. (5) and ( 8 )]. 

We perform a linear least-squares fit of the body observa- 
tions in the detector-network system to the ephemeris posi- 


tions and solve for the registration coefficients (b), see Eq. 
(4). This least- squares fit is performed with the appropriate 
weights, (w rk ,w ck ). 

tek = b 1 + b l( r rpk - (r r *k)) + h(Crpk ~ <<>*»■ 

(4) 

Vek~~ b 4 {^r*k)) b 6 i^rpk (G**))- 

Given the registration parameters of Eq. (4), the tangent 
plane coordinates of the occulted star can be easily calcu- 
lated: 


£* = b I + b 2( r r*k~ ( r r*k )) + h(c r * k -(r r * k ))-bi, 

W*= b 4 + b 5( r r*k" ( r r*k )) + b t( C r*k ” ( C r*k)) = b 4 • 

Now we can calculate the right ascension and declination 
of the occulted star using the transformation from tangent 
plane to celestial coordinates given by Smart (1977): 


tan(a*-a 0 ) = 


cos <%— ? 7 * sin 


tan £* = 


sin <%+ 77 * cos Sq 
cos Sq— rj* sin Sq 


cos^^-ao). 


( 6 ) 


This is the end result of our astrometric reductions. From 
here, the shadow path can be determined by the star’s right 
ascension and declination and the planet’s ephemeris (the 
same one used in the registration above). 


3.3 Calculating the Impact Parameter and Time 

To calculate the impact parameter (or closest approach 
distance between the star and occulting body) we define a 
Cartesian coordinate system (fgh), with its origin at the cen- 
ter of the Earth, with / pointing in the direction of increasing 
right ascension and h pointing to the occulted star (Elliot et 
al. 1993). We construct th efgh coordinates of any object by 
creating its XYZ position from the object’s right ascension, 
declination and distance, d , and then rotating the XYZ coor- 
dinates to fgh. The rotation matrix to convert from XYZ to 
fgh is a function of the right ascension and declination of 
the occulted star (see Elliot et al. 1993). The conversion from 
the (£ 77 ) plane to the (/,g) plane for an object at a distance 
d is f—d^ and g~drj when the same center, (a 0 , 4 )) use( * 
as the origin for each coordinate system. 

To find the geocentric impact parameter, we convert a 
geocentric planet ephemeris (with a light-travel time correc- 
tion) to fgh coordinates. In this coordinate system, we define 
the impact parameter, p, as the minimum distance between 
the center of the occulting body’s shadow and the observer. 
In th efgh system, / O (0 and g 0 (O are the coordinates of the 
observer (Elliot et al. 1993), and f p (t) and g p (t) are the 
coordinates of the center of the occulting body’s shadow as a 
function of time. The impact parameter is given by 

p=Min{\/[/ 0 (f )-/ / ,(/)] 2 + [g 0 (f)-^(0] 2 }- (7) 

The minimization denoted by Eq. (7) gives the impact 
parameter, and the time of the minimum is the predicted 
closest approach time. The latter can be compared with the 
midtime of the occultation light curve as a test ot the astro- 
metric solution. 



3.4 Errors 


Here we calculate the errors in the impact parameter and 
occultation midtime, according to the assumptions of our as- 
trometric reduction procedure. We begin with the errors ex- 
pressed in the tangent plane — derived from Eq. (5), in which 
and 77 * are each expressed as a sum of three terms. There 
is no correlation between the offset term ( b ,) and the slope 
terms ( b 2 and b 3 ) due to the choice of ({r r *k),(c r *k)) as an 
offset for the registration. Also the two slopes terms are as- 
sumed to be uncorrelated because they are in orthogonal di- 
rections. We denote the variance of the sample means 
( r r *k) y {c r *k ) by o- 2 ((r r *k)),cr 2 ((c r *k)) . Taking the appro- 
priate partial derivatives, we find that the variances 
are: 

o 2 (!i*)=o 2 (b i ) + bW((r r * k )) + b]cr((c rn )), 
<r~(v*) = o 1 (t> 4 )+ b 2 5 cr 2 ((r r * k ))+ bl<j 2 ((c r * k )). (8) 

We can gain some understanding of these errors by con- 
sidering a case where the £ axis is perfectly aligned with the 
row axis so that h$ = b 5 = 0. We work only with the equation 
for £, since the 77 result is completely analogous. Using the 
equations for linear-least squares (Clifford 1973), we can 
write an expression for the variance of b i in Eq. ( 8 ) in terms 
of the row weights, w rk , the measurements of the occulting 
body’s position, £ e (t k ) 9 and the scale factor in radians per 
pixel, b 2 . First, we define the ephemeris position in the 
detector-network system: 
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r re (tk)=-j^-+(r r * k ). (9) 

The variance of the sample mean of the star and the 
equivalent quantity for the occulting body (this is equal to 
the variance of the residuals from the ephemeris registration 
expressed in units of pixels instead of radians) are given in 
Eqs. (10) and (11). These are the variances based on the 
scatter; internal errors can be calculated from the weights 
alone, by Eq. (4.31) of Bevington (1992). Note the denomi- 
nators differ because there is one less degree of freedom for 
the occulting body than the star because a line was fit 
through the data, whereas for the star only a mean was de- 
termined. 


° 2 (( r r*k)) = 


t^« r rp k))= 


^i^ifvrW) 2 

(ty- 

(ty-2)2j£ lWjk 


( 10 ) 

(ID 


It is possible that all the frames do not have usable positions 
for both the occulting star and occulted body (for example, if 
the occulting body passes too close to a field star in one 
night’s data). If this is the case, then the number of frames, 
Nj-, for the star and occulting body will not be the same, and 
care is required to implement Eqs. (10) through (12) with the 
proper sums. 

The error in the offset term, b ] is 


J 


o z (b { )=bla 2 ((r rpk )) j -77 


N f^ k = 1 w rk(r rp k-(r r *k)) 2 


j w rk 2 k i , w rk (r rpk -(r r * k )) 2 -[ l w rk (r rpk - (r r * k ))f j ' 


( 12 ) 


If we denote the factor in curly braces by y 2 , then we can 
write (for our ideal case) equations for the variance in both 
£* and tj % : 

) = t>l[ 'f r o 2 ((r rpk )) + <T((r,^))], 

(13) 

^(V*)=b 2 ^.cr((c rpk ))+a 2 ((c r * k ))l 

Equation (13) shows how the errors in the measured occult- 
ing body and occulted star positions translates into errors in 
and 77 * . Both terms in each equation partake of a com- 
mon scale factor (b 2 or b 6 ), and the variance in the sample 
means add quadratically. However, we see that the y factors 
enhance the variance due to the error in the measured body 
positions, increasing as the distance between the mean of the 
body positions is more distant from the star. 

The error in the impact parameter is the error of the star 
position projected onto a line connecting the star at the time 
of closest approach to the center of the occulting body. If 0 is 
the position angle of the star relative to the center of the 
body (measured from North through East) at the time of 
closest approach, then the error in the impact parameter (in 
kilometers) is given by Eq. (14), where d is the distance (in 
kilometers) from the observer to the occulting body. 


aip)=d[o 2 ( £*)sin 2 6+ cr 2 ( 77 * )cos 2 0\ m . ( 1 4) 

Similarly, the timing error is the error component perpen- 
dicular to this, divided by the shadow velocity, v, at the 
midtime of the occultation: 

oito) = - [^( £* )cos 2 <9+ cr 2 (77 % )sin 2 6>] ,/2 . (15) 

4. APPLICATION TO THE TR60 OCCULTATION 

For the Tr60 occultation by Triton, 20 CCD frames of the 
star field were taken on the nights of 1993 July 8-11 with 
the 61 -in. astrometric reflector at the USNO Flagstaff Sta- 
tion. The detector is a Tektronix, thinned, back-side illumi- 
nated, 2048X2048 CCD with an image scale of about 0.325 
arcsec per pixel (field of view equal to 1 1 arcmin). These 
frames were exposed through a broad-band red filter (USNO 
designation A 2 - 1 , with a passband of 6600-8300 A). Typical 
exposure times for the frames were 60 s (longer exposures 
could not be taken because Neptune would saturate). 

Images with known systematic errors were not included in 
the analysis. The data selection process removed centers of 
either field stars, Tr60, or Triton from the analysis; some- 
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Table 1 

Tr60 Astrometric Data, rms Error and Intensity Ratio 


Date # Frames Min. HA Max. HA Tr60 Triton Intensity Ratio 



(Tr60/Triton) 

(hrs.) 

(hrs.) 

Row Error* 

Col Error 8 

Row Error 8 

°i{v)) 

Col Error 8 

4( c n*>) 

(Tr60/Triton) 

1993 July 8 

20/18 

-0.72 

0.62 


0.059 


0.040 

0.5175 ± 0.0006 

1993 July 9 

20/20 

-0.66 

0.90 

0.027 

0.051 


0.059 

0.5267 ± 0.0005 

1993 July 10 

18/2 

-0.81 

0.53 

0.065 

0.076 


0.091 

0.511210.0030 

1993 July 1 1 

16/16 

-0.94 

0.14 

0.074 

0.056 

WESSm 

0.075 

0.4824 ± 0.0009 

a Errors are for an individual frame based on the scatter. Internal error agrees well with error from scatter. Units are pixels and 


the mean focal plane scale is about 0.325 arcsec/pixel. 


times entire frames were removed. Triton centers were re- 
moved either because of high residuals in the ephemeris reg- 
istration (this was the case for two centers) or because Triton 
was merged with a field star (18 of the Triton centers from 
July 10 were not included in the analysis for this reason). 
More detailed modeling of the two sources may provide ac- 
curate astrometry and photometry, but for now these data are 
not included. Two field stars were not included in the analy- 
sis because they were too close to the edge of the CCD (less 
than 20 pixels ~7 arcsec) and showed high residuals in the 
registration. One frame was removed on the basis of photom- 
etry, since the magnitude difference between Tr60 and Triton 
was not consistent with the other frames from that night. 
Frames with thin cirrus as evidenced by inconsistent pho- 
tometry, or field stars not identified on all frames were also 
left out of the reduction. To reduce the effects of refraction 
and telescope flexure, the observations occurred only within 
1 h of the meridian, see Table l. There were 44 field stars 
that comprise the detector network (see Fig. 1), common to 
all the 74 useful frames (Table 1). 


-21:11:40 
-21:13:20 
-21:15:00 

§ -21:16:40 

e 
a 

^ -21:18:20 
- 21 : 20:00 
-21:21:40 
-21:23:20 

19:26:10 19:26:00 19:25:50 19:25:40 19:25:30 19:25:20 

Right Ascension 



Fig. l— Tr60 Field. The star field of the USNO stare frames is shown sche- 
matically. The 44 field stars used to construct the mean-frame coordinate 
system are indicated by filled circles. The occulted star, Tr60, is shown by 
the open circle. The motion of Triton across the field is indicated by the 
diamonds for the first and last observation from each night. The USNO 
frames were all centered to within 3 arcsec. 


4.1 Astrometry 

Before we can apply our astrometric method to the data, 
we must calculate centers for the field stars, Triton and Tr60. 
This is done by putting the flattened frames through a series 
of IRAF scripts that use the functions in the daophot pack- 
age (Stetson 1987) to perform numerical point-spread- 
function (PSF) fitting of the observed stars. First, the routine 
“daofind” is used to identify all stellar images on the frame 
that are within a given range of brightness (neither too faint 
to be in the noise, nor too bright to be in the nonlinear re- 
gime of the detector). This program convolves the data with 
a Gaussian function, and then examines the convolved data 
for local maxima. Approximate centroids, as well as the 
roundness and sharpness of the image are calculated. These 
last two parameters are used to eliminate bright pixels, cos- 
mic rays, and extended objects from the star list. 

Next, aperture photometry of the stars identified by 
“daofind” is performed with “apphot.” An IRAF script se- 
lects stars for the PSF which (i) are neither too bright nor too 
faint, (ii) are not near the edge of the frame, and (iii) do not 
have nearby neighbors. This selection is done on the basis of 
the aperture photometry and “daofind” results. The numeri- 
cal PSF is fit to all identified stars with the daophot func- 
tion “nstar” that uses a least-squares procedure to simulta- 
neously fit stars that are close together. The result is the row 
and column center (r,c), instrumental magnitude, and error 
in the instrumental magnitude for each object. The resulting 
lists are edited to include only: (i) the occulted star (ii) the 
occulting body, and (iii) the detector-network stars. Since the 
errors in the centers were independent of magnitude for the 
network stars, the weighting given by Eq. (1) was used. 

We used the DE21 1 ephemeris to define the motion of the 
Neptune system's center of mass and the NEP 016 model 
(Jacobson et al. 1991) to define the motion of Triton relative 
to the center of mass of the system. Triton’s ephemeris was 
topocentric (i.e., it included effects of geocentric parallax 
and light travel time), but no adjustment was added for the 
apparent displacement due to diurnal aberration, annual ab- 
erration, and refraction. 

One could add the effects of aberration and refraction to 
the ephemeris, but Table 2 shows that the unmodeled contri- 
bution of these effects is negligible. Values of the maximum 
image displacements are given in Table 2. For each effect we 
have selected the frame that would have the maximum image 
displacement and have calculated (i) mean displacements of 
the images for the frame, and (ii) the maximum part of the 
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Table 2 

Maximum Image Displacements in the USNO Frames 


Effect 

Mean for Frame 
(arcsec) 

Not Modeled 
(arcsec) 

Diurnal Aberration 

0.26 

3.4 10- 7 

Annual Aberration 

21. 

2.7 10- 5 

Refraction (monochromatic) 

75. 

6.6 10- 4 

Refraction (chromatic) 

0.015 

8.6 1<H 


effect that would not be modeled by the linear registration of 
Eq. (2). For the effects of diurnal aberration, annual aberra- 
tion, and monochromatic refraction, the appropriate angle 
over which to propagate the unmodeled terms is half of the 
detector field of view. The mean image displacements due to 
chromatic terms in atmospheric refraction do not affect us, 
but the differential displacements due to this effect cause 
unmodeled trends in the positions of network stars that de- 
pends on their color. We have measured the colors for 13 of 
the 44 stars and have calculated the image displacements for 
a hypothetical star with the mean color. The maximum dis- 
placement (determined for the star with the most extreme 
difference in color and for the frame with the greatest hour 
angle) has been entered in Table 2. Although, these image 
displacements would not be removed by the linear registra- 
tion of Eq. (2), most are much smaller than the maximum 
given in the table, and they should be of random magnitude 
over the network. The signature of these errors would be a 
systematic trend of the registration residuals as a function of 
hour angle, but this trend was not apparent. From the entries 
in Table 2 we conclude that the linear approximation for 
these effects would not introduce significant errors into our 
analysis. 

However, if one still wanted to add the effects of aberra- 
tion and refraction to the ephemeris, care is needed to do this 
properly. First one needs to calculate the amount of the dis- 
placement of the Triton center on each frame, then the dis- 
placement added to the ephemeris is the weighted average 
displacements of each individual frame. This is needed to be 
consistent with the Triton centers in the detector network 


coordinate system, which have a net displacement of each 
star by the weighted average of the star’s displacement on 
each frame. 

The residuals from registering the Triton centers to the 
ephemeris are shown in Fig. 2. The figure shows the residual 
of each Triton center as a function of its location in the 
detector network. The residuals have been multiplied by a 
factor of 100. Note that no systematic trends are evident. The 
weighted rms residual is 0.001 arcsec in the £ direction and 
0.003 arcsec in the rj direction. 

The standard error of the Tr60 center and the equivalent 
for the Triton centers (as expressed in the error discussion) 
are listed for each night’s observations in Table 1 . As seen in 
the table, there is a nightly variation of the data quality. 
These numbers also show that the scatter in the Triton posi- 
tion is not inflated due to the presence of Neptune. 

Using the method outlined in Sec. 3, we derived a mini- 
mum distance for the KAO from the center of the shadow of 
359 it 1 33 km (0.017±0.006 arcsec); the center of Triton’s 
shadow was north of the KAO, see Fig. 3. As an independent 
check of the astrometric solution’s accuracy, we can compare 
the observed midtime of the occultation with the closest ap- 
proach time from astrometry. For this event, the astrometric 
closest approach time was 2.2 s before the midtime from the 
light curve. This time difference is less than the formal error 
in the timing (4.1 s) and translates to 60 kilometers (using 
Triton’s shadow velocity of 27.37 kms -1 for the KAO), 
which is less than the formal error on the miss distance. The 
latitudes probed at the half-light level of immersion and em- 
ersion are also listed in Table 3. The full set of latitudes 
probed by this occultation chord are displayed in the right 
panel of Fig. 4. Note that the occultation probes Triton’s 
southern latitudes [those latitudes seen by Voyager (Smith et 
al. 1989)]. The left panel of Fig. 4 shows Triton’s globe and 
the path of the star relative to Triton as seen by the KAO. 

Figure 5 is a plot in the tangent plane of the residuals 
from the registration of the field stars to a secondary astro- 
metric network supplied by A. Klemola (1993, personal 
communication). This network was a re-reduction (using the 
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Fig. 2 — Triton registration residuals. The residuals (multiplied by a factor of 100) from the registration of the Triton centers to the ephemeris are displayed. 
Each residual is indicated by a line whose origin is the Triton position for that frame. 
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Fig. 3 — Triton’s Shadow Path and the KAO. The northern limit, centerline, 
and southern limit of Triton’s shadow are shown as they crossed South 
America. The flight of the Kuiper Airborne Observatory, indicated by the 
bold line, starts at Punta Arenas, Chile and ends in Buenos Aires, Argentina. 
Occultation series images were obtained for 35 min around the occultation 
time on the south-bound leg of the flight. 

ACRS network) of plates taken to search for Neptune candi- 
dates (Klemola and Mink 1991) originally reduced with the 
Perth 70 catalog. The residuals have been multiplied by a 
factor of 100 so that they may be visible in the figure. Note 
that there are no obvious systematic trends in the residuals. 
Since we are only concerned with relative positions, system- 
atic errors in the absolute position (which may not be appar- 
ent in this figure) are not relevant. The standard deviation of 
these residuals are 0.18 arcsec in the £ direction and 0.13 
arcsec in the y direction. 

4.2 Photometry 

The accurate analysis of an occultation light curve relies 
not only on an astrometric calibration but a photometric one 
too. The photometric calibration involves determining the 
relative magnitudes of the occulted star and occulting body 
to define the zero-flux level of the light curve. For the Tr60 
occultation, the relative magnitudes of the objects were de- 
rived from the same USNO data on which the astrometric 
solution is based. There were other sources of photometric 


Tabu: 3 

Tr60 Occultation Circumstances 


Impact Parameter for the KAO(km) 

359 ± 133 

Time Difference 3 (s) 

2.2 ±4.1 

Impact parameter f component (km) 

77 ± 111 

Impact parameter g component (km) 

-350 ± 134 

shadow velocity for the KAO (km/s) 

27.37 

Immersion half-light latitude (deg.) 

-28.9 (+2.6, -2.1) 

Emersion half-light latitude (deg.) 

3.6 (+4.5, -5.2) 


a Light curve occultation midtime minus astrometric closest 
approach time 


data, but due to nonlinearities in other detectors, these data 
are the best source for the photometric calibration. 

The flux ratios of Tr60 to Triton on a nightly basis are 
listed in Table 1 and are displayed in Fig. 6. There is vari- 
ability in the relative intensities, either from variability in 
Triton as a function of rotational phase or from the variabil- 
ity of Tr60. To get the intensity ratio at the occultation time, 
we linearly interpolated using the last three nights of data (as 
defined by the midtime of the light curve). This photometric 
determination of the relative intensities involves a different 
instrument and filter than the occultation observations. The 
effective wavelength of each system (KAO and USNO) was 
determined from the quantum efficiency of each detector and 
the filter (no filter was used for the KAO observations). 

The effect of refraction on the effective wavelength was 
ignored. The constant of refraction changes by less than \% 
over the wavelength of the filter used at the USNO. At the 
altitude of the KAO, the atmospheric pressure is lower than 
sea level by about a factor of 5, and the change in refraction 
over the wavelengths of the instrument (the NASA/Ames 
Portable CCD) sensitivity is only 1%. R and / photometry of 
Tr60 (/? = 13.39±0.02; /= 13.82±0.02; D. M. Chamberlain, 
personal communication) and Triton (Buratti et al. 1994) 
provided the slope of the intensity ratio with wavelength 
from the effective wavelength of the USNO observations 
(0.65 yLtm) to the effective wavelength of the KAO observa- 
tions (0.67 /^m), see Table 4. The R and / magnitudes for 
Tr60 are approximately Kron-Cousins (M. Buie, personal 
communication). 

5. DISCUSSION 

The method outlined in this paper to determine the closest 
approach distance between a star and a body passing close by 
has been applied to the 1993 Triton occultation, and we find 
consistency between the closest approach time predicted by 
the astrometry and that derived from the occultation light 
curve. To the level of accuracy presented, we find no system- 
atic problems with either the analysis or the data. Much of 
this is due to the high astrometric quality of the telescope and 
detector, as well as the practice of observing only when the 
objects are near transit. 

Now we consider the sources of formal error. Note that 
choosing the occulted star for the origin of the detector- 
network coordinate system simplifies the error analysis, so 
that only the errors in the offset terms of the registration and 
the scatter of the star center [the two parts of Eq. (13)] con- 
tribute. For this dataset, we are limited by the error in the 
registration coefficients (af6|)= 0.005 arcsec and 
oib 4 )=0.006 arcsec), not the scatter in the Tr60 center 
(0.001 arcsec in £ and 0.002 arcsec in y). Therefore, more 
observations of the star would not significantly improve the 
solution. Improving the precision of the solution could have 
been accomplished by (i) more observations of Triton (to 
reduce the registration coefficient error), which cannot be 
increased arbitrarily because the occulting body only spends 
a finite amount of time in the same field with the occulted 
star, and (ii) arranging the observations so that the mean 
Triton position were close to the Tr60 position. If this were 
true for the present data set, then the y factor in Eq. (13) 
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Time Past 22:54:15 UT (Sec) 


Fig. 4 — The region of Triton probed by this stellar occultation. The left panel shows Triton as seen by the KAO. Triton’s south pole is visible. The path of the 
star relative to Triton is indicated at 5-s intervals. The right panel shows the latitude of Triton that was probed as a function of time. The dashed lines indicate 
the ± 1 -sigma limits on the latitudes probed. Note that the time axis is reversed to be consistent with the path of the star in the left panel. The midtime of the 
occultation, along with the half-light immersion and emersion times from the light curve are indicated. The occultation probed the southern hemisphere during 
the midtime of the occultation. 


would have been unity and the error in the minimum sepa- 
ration would have been reduced from 0.006 to 0.003 arcsec. 
For this analysis, the effective mean Triton position was fur- 
ther from the Tr60 position than a simple (unweighted) av- 
erage would indicate, because the first two night’s data (both 
recorded before the event) were the highest quality. 

Outer solar system bodies such as Triton and Pluto are on 
the same 1 1 arcmin field as a given candidate star for about 
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Fig. 5 — The residuals from registering the weighted mean frame to a sec- 
ondary astrometric network. The residuals have been multiplied by a factor 
of 100. Each of the standard stars are indicated in the figure by their (£77) 
components, and their residual components are indicated by the arrow. 


one week at opposition. However, less remote bodies gener- 
ally will move through the field more quickly, decreasing the 
window of opportunity for obtaining this critical calibration 
data. For the most accurate results, the observations should 
be taken both before and after the occultation, so that we can 
interpolate to derive the closest approach distance (rather 
than extrapolate) and to minimize the y factor in Eq. (13). 



Fig. 6 — The relative intensity of Tr60 and Triton as a function of time. The 
filled circles indicate the intensity ratio for each night’s observations with an 
error bar given by the standard error of the sample mean. There were only 
two good Triton observations on the third night, hence the error bar is larger. 
Also shown in the figure is the best fit line to the observations from the last 
three nights. The open circle represents the adopted intensity ratio measured 
with the USNO system for the time of the occultation. 
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Table 4 

Tr60 Photometry 


Dataset 

effective X 

Intensity Ratio 


(Hm) 

(Tr60/Triton) 

USNO 

0.67 

0.4910 ±0.0007 

KAO 

0.65 

0.4895 ± 0.0012 a 


a derived from the USNO intensity ratio and 
the effective wavelengths 


This method provides the following two advantages in 
comparison to the standard astrometric reduction technique 
which is commonly used for occultation predictions (Wasser- 
man et al. 1979). First, no astrometric network is needed. 
With CCD data, the field is generally small enough that pri- 
mary astrometric networks do not have a sufficient number 
of standard stars. So eliminating the need for the dense as- 
trometric network saves significant time and effort. The sec- 
ond advantage is that this method is immune to scale and 
rotation errors in an astrometric network that result from the 
random errors in the star positions and any systematic differ- 
ences between the stellar reference frame and that of the 
ephemerides. We do not expect to find significant scale errors 
in the planetary motion (which sets the scale) as defined by 
ephemerides, because this error would accumulate over time 
and an offset of the body from its ephemeris would be no- 
ticed. 

This method has the disadvantage that it cannot be used 
for all occultations. If the occulting body has a linear path 
across the detector, then the motion of the body can only set 
the scale in the direction of motion, and an astrometric net- 
work is needed to define the perpendicular scale. This “com- 
promise” method was used in the 1994 March 9 occultation 
by Chiron (Elliot et al. 1995). 

Over time, this ephemeris method can be employed to 
determine a “catalog” of star positions in the reference 
frame of planetary ephemerides. This catalog could then be 
used to investigate differences between the reference system 
of ephemerides and stellar reference frames. Toward this 
end, we present three positions for Tr60 (all in J2000) based 
on three different astrometric reductions, see Table 5. These 
positions come from the ephemeris method and from astro- 
metric reductions of plates (Klemola 1993, personal commu- 
nication) using two separate primary astrometric networks: 
the ACRS and Perth 70. The coordinates of the detector- 
network stars are not given here because the right ascension 

Table 5 

Tr60 Positions (J2000 FK5) 


Reference System RA Dec 

Triton ephemeris 19 25 37.227 -21 17 32.66 

ACRS network 19 25 37.256 -21 17 32.70 

Perth 70 network 19 25 37.205 -21 17 32.73 


and declination of these stars were not determined as accu- 
rately as Tr60 since Triton did not pass close to these other 
stars. Many of the stars would be outside of Triton’s path, 
and the image scale would have to be extrapolated to their 
locations, which introduces errors. 

For the most precise determination of an occultation 
chord location, it is best to have well spaced, multiple 
chords. However, with occultations by small bodies it is not 
always possible to have multiple observers in the shadow 
path (due to the uncertainties of the prediction, even a week 
in advance). Useful information can be derived from these 
single-chord occultation observations with the ephemeris 
method presented here, if there are quality astrometric data 
of the occulted star and occulting body both before and after 
the occultation. 
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